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ABSTRACT 

Recent spectroscopic models of active galactic nuclei (AGN) have indicated that the recommended electron- 
ion recombination rate coefficients for iron ions with partially filled M-shells are incorrect in the temperature 
range where these ions form in photoionized plasmas. We have investigated this experimentally for Fe XIV 
forming Fe XIII. The recombination rate coefficient was measured employing the electron-ion merged beams 
method at the Heidelberg heavy-ion storage-ring TSR. The measured energy range of 0-260 eV encompassed 
all dielectronic recombination (DR) ls 2 2s 2 2p 3131' 31" nl'" resonances associated with the 3p\/2 — ► 3/?3/2, 
3s — > 3p, 3p — > 3d and 3s — > 3d core excitations within the M-shell of the Fe XIV (ls 2 2s 2 2p 6 3s 2 3p) parent ion. 
This range also includes the Is 2 2s 2 2p 6 31 31' 41" nl'" resonances associated with 3s — » 41" and 3p — > 41" core 
excitations. We find that in the temperature range 2-14 eV, where Fe XIV is expected to form in a photoionized 
plasma, the Fe XIV recombination rate coefficient is orders of magnitude larger than previously calculated 
values. 

Subject headings: atomic data — atomic processes — plasmas — galaxies: active — galaxies: nuclei — 
X-rays: galaxies 



1. INTRODUCTION 

Recent spectroscopic XMM-Newton and Chandra X-ray 
observations of active galactic nuclei (AGN) have detected 
a new absorp tion feature around 15-17 A JSako e t alJ|200l|: 
Pounds et al] l200lt iKalpi et alJ l200l IBehar et alJ 2003[ 
Steenb Tugge et al j|2003t iKasni et al.H2004UGallo et alJl2004 
Pounds et al. 2004; Krongold et al. 2005). This has been iden- 



tified as an unresolved transition array (UTA) due mainly to 
2p — > 3d inner shell absorption in moderately charged iron 
ions with an open M-shell (Fe I-Fe XVI). On the basis of 
atomic structure calculation s and photoabsorption modeling, 
IBehar. Sako. & Kahnl (120011) pointed out that the shape of the 
UTA features can be used for diag nostics of the AGN ab- 
sorber. However, Netzer et alj fc003l) noted a disagreement 
between the predicted and observed shape of this feature. As 
a possible cause for this discrepancy they suggested an under- 
estimation of the low temperature dielectronic recombination 
(DR) rate coefficients for iron M-shell ions. These rate coef- 
ficients determine the charge state balance of the iron M-shell 
ions in a photoionized plasma and, consequently, the shape of 
the UTA. 

Reliable low temperature DR rate coefficients of iron M- 
shell ions are not availabl e in the literature. The widely used 
compilation of Arnaud & Raymond ( 1992) is largely based 
on theoretical work by Jacobs et al. ( 1977) and Hahn ( 1989). 
The purpose of this early theoretical work was to produce DR 
data for modeling coronal equilibrium. However ions form in 
coronal equilibrium at temperatures about an order of magni- 
tude higher than those where they form in photoionized gas 
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llKallman & Bautistall2001l) . Therefore, it is questionable to 
use these theoretical DR data for photoabsorption modeling. 
Benchmarking by experiment is highly desirable. 

The only DR measurements for M-shell ions of third- and 
fourth-r ow elements available up to now are th ose for Na-like 
Fe XVI JLinkemann et al.l l 1 995t iMiillerll 1 999t) : Ar-like Ti V 
and Sc I^S chippers "alll99al27)02lh and Mg-like Ni XVIII 
(Foale et al. 2003). Although being of minor astrophysical 
relevance, the Sc IV and Ti V measurements illustrate that 
the low-energy recombination spectra of M-shell ions can be 
dominated by strong resonances associated with 3p — ► 3d core 
excitations. 

This work presents the experimental measurement of ra- 
diative recombination (RR) + DR for Fe XIV forming 
Fe XIII. In the temperature range where Fe XIV is ex- 
pected to exist in a photoionized plasma, our experimentally- 
derived RR + DR rate coefficient is orders of magnitude 
larger than the sum of the r ecommended RR value from 
Wood s. Shull. & SarazTnl ill 9811) and the P R value from the 
compilation of Arnaud & Ravmondj lll992l) . As DR is much 
larger than RR at these temperatures, this discrepancy is 
clearly due to errors in the recommended DR data. A sim- 
ilar discrepancy can be expected for the other iron M-shell 
ions. This has important consequences for the modeling of 
AGN spectra. 

2. EXPERIMENTAL TECHNIQUE 

The experiment was performed at the heavy-ion test stor- 
age ring TSR of the Max-Planck-Institut fiir Kernphysik in 
Heidelberg, Germany (MPI-K). Measurements employed the 
well established procedures for studying electron-ion recom- 
bination. Details of experimental and data r eduction proce- 
dures as well as f urther references ar e given in Schippe rs et alJ 
j200lll2004l) and lSavin et afllEool . 

A beam of 56 Fe XIV was provided by the MPI-K acceler- 
ator facility. After acceleration to an energy of 4.2 MeV/u 
the ion beam was injected into the storage ring where it was 
collinearly overlapped with a magnetically guided electron 
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FIG . 1 . — Measured Fe XIV to Fe XIII merged-beams electron-ion recom- 
bination rate coefficient (MBRRC) in the energy range associated with DR 
resonances via 3/^/2 — > 3* ~ * 3p and 3p — > 3d core excitations. The 

vertical bars denote DR resonance positions as expected on the basis of the 
hydrogenic Rydberg formula. Only the (3s 2 3d ^D^j^nl'" series of Rydberg 
resonances can unambiguously be identified. Note that the resonances below 
2.3 eV exceed all other resonances in height by an order of magnitude. 
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FIG. 2. — Measured Fe XIV to Fe XIII merged-beams recombination rate 
coefficient (MBRRC) at energies below 10 eV (filled circles). The full curve 
enclosing the white area is the result of a fit comprising of 39 DR resonances 
and the contribution from RR (short-dashed curve). The dashed curves show 
the fits to the first 6 of the 39 DR resonances. The excess rate coefficient 
Aa (see text) contributes to the measured signal at energies below Ri 1 meV. 
Here the assumed enhancement factor is (An/an) + 1 = 1.7 



beam. In the overlap region the ions can recombine with 
the electrons via RR and DR. Recombined ions were sepa- 
rated from the ion beam in the first dipole magnet after the 
electron-ion interaction region and counted using a single- 
particle scintillation-counter with a nearly 100% detection ef- 
ficiency. The systematic experimental uncertainty of the mea- 
sured recombination rate coefficient is estimated to be ±15% 
(Lamp ertet alJl9 96). This uncertainty stems mostly from the 
ion current measurement. Adding further uncertainties dis- 
cussed in Sec. [3]increases this to 18%. 

The experimental electron energy distribution is best de- 
scribed as a flattened Maxwellian distribution which is char- 
acterized by the longitud inal and transver se temperatures 7]i 
and T± (see, e.g., Pastu szka et aljfl996h . The experimen- 
tal energy spread AE = [(ln(2) k B T ± ) 2 + 161n(2)£fe B 7||] 1/2 de- 
pends on the electron-ion collision energy. Here Icq is the 
Boltzmann constant. In the present experiment the temper- 
atures were k^T± ps 12 meV and ksTu m 0.09 meV and con- 
sequently AE = 8.3 meV for E < 0.2 eV, AE = 33 meV at 
E = 1 eV and AE = 245 meV at E = 60 eV. In the following 
the experimental data are presented as a merged-beams re- 
combination rate coefficients (MBRRC) (crv), i.e., the recom- 
bination cross section times the relative velocity convolved 
with the experimental electron energy distribution function. 

3. RESULTS AND DISCUSSION 

Figure \l\ shows the measured Fe XIV MBRRC in the en- 
ergy range of 0-60 eV. This includes all resonances associ- 
ated with 3s 2 3p — > 3a ,2 3Z" and 3s 2 3p — > 3s3p3l" core exci- 
tations. The experimental MBRRC increases strongly at en- 
ergies below 2.5 eV. Figure [2] shows this energy region in 
more detail. It is dominated by strong DR resonances that 
presumably represent levels of the type (3s3p 2 2 Pi/2, 3/2)7/"'. 
Furthermore the high-/ limit of the (3s 3p 2 2 Pi/2. 3/2)7/"' res- 
onances a nd the series limit o f the 3s 2 3p3/2«/'" series at 
2.337 eV (Behring et al. 1976) are quasi degenerate within 
wO.l eV. DR resonances associated with a 2p — * Pi/2 
fine structure excitation were found to be important for low- 



temperature DR of Fe XVIII forming Fe XVII (ISavin et alJ 
I1997LI1999I) . The measured Fe XVIII MBRRC exhibits a reg- 
ular Rydberg series of (2p 5 2 Pi/2)«/"' resonances. In the 
present Fe XIV case such a regular resonance pattern from 
the 3/?3/2«/"' series is not observed. The question remains 
to what extend the unusually strong rise in DR resonance 
strength below 2.5 eV is caused by the near degeneracy of 
(3s3p 2 2 Pi/2. 3/2)7/"' and 3s 2 3/?3/2«/'" resonances. Accu- 
rate theoretical calculations taking into account these features 
and their mutual coupling through quantum mechanical mix- 
ing appear highly desirable in order to understand the detailed 
electronic dynamics leading to the high resonant rates in this 
particular energy range. 

Up to an energy of w 50 eV the resonance pattern is 
irregular. Assignment of the observed features could not 
be achieved with the exception of the (3s3p 2 2 S)n/"' and 
(3s3p 2 2 P)nl"' series limits around 45 and 48 eV, respec- 
tively. At higher energies the (3s 2 3d 2 Dj,/2)nl"' Rydberg 
series is discernable converging to its limit at 58.6722 eV 
(Behrinsetal. 1976). This feature was used to calibrate the 
experimental energy scale. This was achieved by multiplying 
the experimental energy scale by a factor that deviates from 
unity by less than 1%. 

Finally, Fig. [3] shows the measured Fe XIV MBRRC in the 
range of 60-260 eV. This range includes resonances associ- 
ated with 3s — > 41" and 3p — > 41" (i.e., AN = 1) core excita- 
tions. Here N is the principal quantum number of the transi- 
tioning core electron. The strongest features in this spectral 
range are the series limits of the 3p — > 41" nl'" resonance se- 
ries. These occur at about 177.9 eV (/" =0), 195.0 eV (/" = 1), 
210.4 eV (/" = 2), and 22 1.7 eV (/" = 3) accordi ng the NIST 
atomic spectra data-base (Ralchen ko et alJl20 05). Identifica- 
tion of the observed peak features with individual Rydberg 
series is challenging except for the series limits. 

The plasma recombination rate coefficient (PRRC) is de- 
rived by convolving the measured MBRRC with a Maxwell- 
Boltz mann electron energy distribution. As detailed by 
Scrap pers et alJ d2001l [2004), there are three issues that re- 
quire special consideration: the experimental energy spread, 
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TABLE 1 
Parameters for the fit of Eq.Q 
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3. — Measured Fe XIV to Fe XIII merged-beams recombination 



rate coefficient in the electron-ion collision-energy region of the 3.v 2 4/"n/'" 
and 3s 3p4l" nl'" resonances attached to the 3p — » 4/" and 3s — ► 41" core 
excitations, respectively. 



the recombination rate enhancement at low energies, and field 
ionization of high Rydberg states in the storage-ring bending 
magnets. 

The experimental energy spread AE influences the out- 
come of the convolution for resonances with resonance en- 
ergies E < AE. This can be circumvented by extracting the 
DR cross section from the measured MBRRC at low ener- 
gies. In Fig.|2]this was achieved by fitting 39 DR resonance 
line-shapes to the measured spectrum. It should be noted that 
most probably none of the 39 DR line-shapes correspond to 
a single resonance but each rather comprises a blend of res- 
onances. Additionally the fit includes a contribution due to 
RR which was calculated using a semi-e mpirical hydrogenic 
formula (Eq. 12 offSchippers et al.l2001l) . 

An enhanced MBRRC is consistently observed in merged 
electron-ion beam experiments at very low energies. For 
E < kftTu the measured MBRRC excee ds the theoretical ex - 
pectation by factors of typically 2-3 (Gwinner et al. 2000). 
This excess rate coefficient (labeled Act in Fig.|2j is an arti- 
fact of the merged-beams technique and hence a mean excess 
rate coefficient was subtracted prior to the derivation of the 
PRRC. 

Field ionization of the loosely bound high Rydberg electron 
in the recombined ions can result from the motional electric 
fields they experience inside the storage ring bending magnets 
(Schippers et al. 2001). In the present experiment only RR 
and DR involving capture into Rydberg levels with quantum 
numbers less than 55 contribute to the MBRRC. 

Similar to the approach of Schippers et al. (2001, 120041) and 
iFoele et aTl (120051) . the missing DR resonance strength up to 
"max = 1000 was estimated from a theoretical calculation us- 
ing the AUTOSTRUCTURE code dBadnelll[i 986). Although 
the AUTOSTRUCTURE code fails to reproduce the irregu- 
lar Fe XIV resonance structure below 50 eV, it reproduces 
the more regular resonance structures of high-n Rydberg reso- 
nances close to the various series limits reasonably well when 
slight 'manual adjustments' are made to the autoionization 
transition rates. 

The unmeasured DR contribution due to n > 55 was esti- 
mated to contribute (15 ±7)% to the experimentally-derived 
Fe XIV PRRC for plasma temperatures k B T e >9 eV. The 7% 
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uncertainty was estimated from the adjustment of the transi- 
tion rates. The same field ionization cut off was taken into 
account in the RR contribution to the low energy resonance 
fit in Fig. The unmeasured fraction of the RR-PRRC, due 
to RR into 3s 2 3pnl (n > 55) states, was calculated using a 
hydrogenic formula (Eq. 13 of Schippers et al. 2001, Fig. @}. 
The RR and DR contributions from n = 55-1000 were subse- 
quently added to the PRRC. The contribution of AN = 1 DR 
with n > 55 is insignificant. 

For convenient use in astrophysical modeling codes the to- 
tal Fe XIV to Fe XIII PRRC was fitted using 



aplasma^e) = T e 3/2 V] Ci&X^-Ei / k B T e ) 



i=l 



(1) 



The resulting fitting parameters c, and £, are given in Table 
[l] The fit deviates by less than 1 % from the experimentally- 
derived result in the energy range 70 meV- 10 keV. 

Uncertainties in the non-resonant portion of the background 
that was not due to RR contributed an estimated 8% error to 
the measurement. Uncertainties in the exact enhancement fac- 
tor near eV contributed a 2.1% uncertainty at 0.07 eV, and 
1.4% at 0.1 eV, and less than 1% above 0.14 eV. Adding 
these, the transition rate uncertainty, and the 15% system- 
atic uncertainty in quadrature gives a total experimental un- 
certainty of ±18% in the energy range 70 meV- 10 keV. 

4. ASTROPHYSICAL IMPLICATIONS 

The approximate temperature ranges where Fe XIV forms in 
photoionized and in co llisionally ionized plasmas can be ob- 
tained from the work of Kallman & B autista ( 200 1 ) . For pho- 
toionized plasmas they find that the fractional Fe XIV abun- 
dance peaks at a tem perature of 4.5 eV. From a differe nt pho- 
toionization model, Kraemer. Ferland. & Gabel (2004) obtain 
a value of 7.8 eV. The 'photoionized zone' may be de- 
fined as the temperature range where the fractional abun- 
dance of a given ion exceeds 10% of its peak value. For 
Fe XIV this corresponds to a temperature range of 2-14 eV 
( Kallman & Bautista 2001). Usin g the same criterion and 
again the results of Kallman & Bautista] ( 1200 ll) . for coronal 
equilibrium the Fe XIV 'collisionally ionized zone' is esti- 
mated to extend over a temperature range of 170-580 eV. It 
should be kept in mind that these temperature ranges are only 
indicative. In particular, they depend on the accuracy of the 
underlyingatomic data base. 

In Fig. |3] we compare our expe rimentally-derived P RRC 
with the RR rate coefficient of IWoods. Shull. & Sa razin 
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Electron temperature (eV) 

FIG. 4. — Experimentally-derived Fe XIV to Fe XIII recombination 
rate coefficient in a plasma (PRRC, thick full line) comprising AN = DR 
(Fig.QJ, AN = 1 DR (Fig.|3j, RR, and the theoretical estimate for the un- 
measured contributions of states with n > 55 for RR and AN = DR. The 
error bars denote the ±18% experimental uncertainty in the absolute rate co- 
efficient. The experimental results without the RR and DR extrapolation are 
shown by the dash-dot-dotted line. Also shown is the theoretical DR rate 
coefficient of Arnaud & Raymond 1 1992, thick dashed line, labeled AR), its 
deliberate modification by Netzei 1 2004 , dotted line), and the RR calculation 
bv lWoods. Shull. & Sarazinl <198ll mln dotted line, labeled WSS). The sta r 
at 4 eV (labeled KFG) is the estimate of lKraemer. Ferland. & Gabel l200l) . 
The contribution from the low energy DR resonances between and 2.5 eV 
is shown as the thin dashed line. T he curves labeled RR were calculated us- 
ing a hydrogenic formula (Eq. 13 of Schippers et al. 2001) with ;i max = 1000 
(thin full line) and with n raax = 55 (dash-dotted line). The temperature ranges 
where Fe XIV is expected to peak in abundance in photoionized and colli- 
sionally ionized plasmas are highlighted. 



ill 98 11) and th e recommended DR-PRRC of 

Arnaud & Rav monJ d!992|). This l a tter is based on a 
theoretical calculation by iJacobs et alJ i ll 9771) and includes 
DR associated with 2p — > 3d inner shell transitions. It was 
calculated for a temperature range of k%T = 30-9000 eV. 
In the collisionally ionized zone the experimentally-derived 
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PRRC is lower than the theoretical result by 21% at 
k^T = 580 eV. This discrepancy may partly be attributed to 
the fact that DR resonances associated with 2p — + 3d and 
higher core excitations were not measure d. But given t h e poor 
agreement between the calculations of IJacobs et alJ (119771) 
and P R measurements for other iron ions (e.g., Savin et al. 
2002), this close agreement is probably fortuitous. 

In the photoionized zone the experimentally-derived PRRC 
is decisively determined by the resonances occuring at 
electron-ion collision energies below 2.5 eV. The RR con- 
tribution is insignificant relative to DR. Our derived PRRC 
is several orders of magnitudes larger than the DR data from 
the widely used compilation of Arnaud & Raym ond! lll992l) . 
They are still ab out an order o f magnitud e larger than t he DR- 
PRRC used bv lNetzerl i l2004l) and lKraemer. Ferland. & Gabell 
(2004) who deliberately assumed sets of DR-PRRC for the 
Fe M-shell ions in order to achieve a better agreement of their 
plasma modeling results with the observed shape of the Fe 
2p — > 3d UTA (cf. introduction). 

The present result shows that the previously available the- 
oretical DR-PRRC for Fe XIV forming Fe XIII are much too 
low, as Netzer et al. (2003) had already suspected. Other stor- 
age ring measurements show similar deviations from pub- 
lished recommended low temperature DR-PRRC for M-shell 
ions ( iLinkemann et alJl995tMulleill999t Fogle et alJ2003l) . 
We are now in the process of carrying out DR measurements 
for other M-shell iron ions. As they become available we rec- 
ommend that these experimentally-derived PRRC be incorpo- 
rated into future models of AGN spectra in order to arrive at 
more reliable results. 
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